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Abstract 

We present a simple method to calculate certain sums of the eigen- 
values of the volume operator in loop quantum gravity. We derive the 
asymptotic distribution of the eigenvalues in the classical limit of very 
large spins which turns out to be of a very simple form. The results can 
be useful for example in the statistical approach to quantum gravity. 

1 Introduction 

Volume operators were introduced by Rovelli and Smolin £Q and by Ashtekar 
and Lewandowski [2] (the latter will be used in the current paper). The 
volume operator is much less understood than the area operator, one of the 
reasons being that it is not diagonal in any reasonable basis and there are 
ambiguities in defining its action on a given graph (see |3j and references 
therein; clarification of the issue is given in The matrix elements of the 
volume operator were obtained in [3], [U| (in the context of cosmology [7j); 
however, the formulae are very complicated and it is extremely difficult to 
get analytic expressions for the eigenvalues of the volume operator, except 
for very small values of the momenta (for some bounds on the eigenvalues 
see Therefore it seems worthwhile to try to find some alternative way 
of describing the eigenvalues. In this paper we show that relatively simple 
formulae can be obtained for the sums of squares (or higher even powers) of 
the eigenvalues, with the most important "quantum numbers" in the vertex 
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fixed: j, jx, j 2 , js, where j is the total spin and ji, j 2 , J3 are spins of the 
incoming legs. We also calculate the sums of higher powers (4, 6 and 8) of 
the eigenvalues with fixed ji,j 2 , J3- At the end of the paper we derive (in the 
form of an integral) the distribution of the eigenvalues for fixed 31,32,33,3 
in the classical limit, when all the values of spins tend to infinity, and show 
that without fixing j the integral can be performed; the result turns out to 
be surprisingly simple: p(X) = arccos(\/ (j 1 j 2 j 3 ) / (2j 1 j 2 j 3 ) for |A| < jij 2 h 
and p(X) = otherwise. 



2 Volume operator 

The volume operator V is defined as 

V = \M (1) 
where q is analogous to the determinant of the metric 

q= e ijkJ} Jj Jr- ( 2 ) 

K.J<K 

Here, /, J, K run through 1, 2, N (N > 3) and k are SU{2) indices and 
we omit an overall constant proportional to the Planck volume. Eigenvalue 
of V describes a contribution to the volume from a given vertex with N + 1 
legs with "spin conservation'. We will concentrate on the most important 
case iV = 3 (so called 4-valent vertex) and later we will comment on the 
application of the method to higher N. 

The operator q has several important properties: 

• It commutes with all J 2 j 

[9,Jj] = 0. 1 = 1,2,3 (3) 
so we can treat in the following all Jj = + 1) as fixed. 

• It commutes with all components of the total spin J 

[g,J]=0. J:=J! + J 2 + J 3 (4) 

so its eigenvalues can be further indexed by the eigenvalues of J 2 , de- 
noted by j(j + 1). Note that because of (@J the eigenvalues of q 2 are 
the same for all m in the (j,m) multiplet. By this procedure, we are 
automatically in the gauge invariant case (where J 4 := — J satisfies 

Jl + J 2 + J 3 + J4 = 0). 
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operator q treats all four spins Ji, J 2 , J3, J4 satisfying J x + J 2 + J3+ J4 = 
in the same way (modulo sign) since for example 

Q = £ijkJ\ J2 J3 = £ijkJ\ J2 J\ (5) 



• The operator q therefore suggests a natural way to number states in a 
vertex: we specify ji, 32, h, j, Tn of incoming legs and the eigenvalue 
of q (the numbering is usually done in a "non-democratic" way: spins 
ji and j 2 are composed to spin j, spins j 3 and j 4 are composed also to 
spin j and the two spins j are composed to a singlet) 

Since the operator q is self adjoint and (as can be easily seen from anti- 
symmetry of €ijk and cyclic property of trace) trace of any odd power of q 
vanishes, therefore q has only real eigenvalues that are either or come in 
pairs (Aj, — Aj). We assume everywhere that all ji are bigger than since if 
any of jj is equal to then all the eigenvalues of the volume operator are 
also equal to 0. 

As was shown in jH] , jB] operator q has very complicated matrix elements as 
a function of 3,31,32,32, (the matrix elements were obtained from the formula 
q = i[Ji • J2, J2 • J3] and projection on a state with given j). 

Note that there exists a convenient representation of the operator q in 
the space of polynomials of degree (2j 1 + 1) in x, (2j 2 + 1) in y and (2j 3 + 1) 
in z where 

\j u mi ) = ^(2 Jl )\/( Jl - mi )\/( Jl+mi )\ x n+^ ( 6 ) 
and the generators 

Ji = -x 2 d x + 2j x x, Jf = d x , Jj 5 = xd x - ji (7) 
and analogously for two other sets of generators. Then q has the form: 

q = ^[-(x-y)(y-z)(z-x)d x d y d z +j 1 (y-z)(y + z-2x)dyd z 

+j 2 (z -x)(z + x- 2y)d z d x + j 3 (x - y){x + y - 2z)d x d y 
+2j 1 j 2 (x - y)d z + 2j 2 h(y - z)d x + 2 ]3]l {z - x)d y ] (8) 

As we can see explicitly, the q operator conserves the total degree of the 
state (i.e. m = mi + m 2 + m 3 ) and indeed has complicated matrix elements 
especially after projection on eigenstates of J 2 . 
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3 The square of the determinant operator 

We consider the second power of the operator q: 

ffi — p P T k \ 2 

H — \ t ijk J l J 2 J 3 J ■ 

Using the well known formula 

tijk^lmn ^il^jm^kn ^jl^im^kn ^kl^im^jn 

^il^km^jn ^jl^km^in ^kl^jm^in 

we obtain 

q 2 = J^JgJ^ — J x J|(J2 ■ Js)^3 — ' J2)^2 ~~ ^2^1(^1 ' ^3)^3 



(9) 



(10) 



+Jl(j 1 • j 2 )(j 2 • j 3 )j^ + J[J 3 2 {h ■ j 3 )4j 3 3 . 

We want to trace the operator with the states 

|(j'i,mi) ® (j 2 ,m 2 ) ® {h,mz)) 



(11) 



(12) 



(often written as |mi,m2, 777.3) for short), not yet imposing any restriction 
on j and m = rri! + 777 2 + 777.3. The calculation uses well known formulae 
(straightforward from (jUJ) and (|7j)): 



Jf\ji, ™i) 
Jj\ji, rnj) 
Jflji, ™<i) 
The result reads 



jl(jl + 1) - 777/(777/ + 1) 777/ + 1) 



3l(jl + 1) - 777/(777/ - 1) |j/, 777/ - 1) 
777/ |j/, 777/) 



(13) 



(mi, 777 2 , 777 3 | q 2 | 777i, 777 2 , 777 3 ) 
1 



-777 



- 1 

1 2 
+ -777 9 

2 2 

1 2 
+ -777o 

2 3 

Summing over all mi, m 2 , 777.3, we have the result 
Tr g 

y 
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(J2O2 + 1) - ^2) (j3(j3 + 1) - ml) - m 2 m 3 

(j'l(jl + 1) - m 2 ) (j 3 (j 3 + 1) - 777 2 ) - 777x7773 
(ilO'l + 1) - ^1) (j2(j2 + 1) - 777 2 ) - 777i777 2 



2 f y n(, n + l)(2ji + l )./,(.,, + l)(2j 2 + l)j 3 (j 3 + \ )(■>,, + 1). 



(14) 



(15) 



Since the total number of states is equal to (2ji + 1) (2j 2 + 1) (2j 3 + 1) the 
average eigenvalue of q 2 is equal to 

(f) = ihUl + 1)J2(J2 + l)h(j3 + I)- (16) 

In the last section we will derive the general formula for (q 2N ) in the limit of 
large j's. 



4 Sum of the eigenvalues with fixed j 

In this section we calculate Tr q 2 with an additional condition 

m i + m 2 + m 3 — m, (17) 

where to be definite we take m < 0. We sum in this way many different j 
multiplets with j > — m; it is important to notice that because of (jSJ) the 
eigenvalues of q 2 are the same for all (j, m) components of a given j multiplet. 

As one can easily convince oneself it is a non-trivial problem to impose 
the condition (j!7|) directly Therefore we use an indirect method that will 
partly do the bookkeeping for us: we perform the sum with the additional 
factor 

f2(w,j 1 ,j 2 Js) : Tr (q 2 u->--> (18) 

and at the end we will get to the final result by picking up the coefficient in 
front of w m . 

We use the formula 

£ m p lW ^ = (wd w Y W . (19) 

mi=-ji W 

The computation of f2(w, ji, 32, J3) from (fTIjl gives 

2 

(1 - w)f 2 {w,j 1 ,j 2 ,j3) = j-, vg w -n-T2-ra+i {j 1 j 2 j 3 (j 1 + j 2 + j 3 ) 

(1 — ifj 8 

-w [jiMUi + j 2 )(2j 3 + 1) + 2j 3 - 1) + (1 <- 3) + (2 ^ 3)] 
2 

w 3 [(ii + l)(j 2 + l)(0i + j*)(2j 3 + 1) + 2j 3 - 1) + (1 ^ 3) + (2 ^ 3)] 



+w 2 



(Ji + l)0'i + 5O2 + h)(2j! - 1) + 6JU2J3) + (1 » 2) + (1 ^ 3) 
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+ W \ 3l + l)(j 2 + l)(j 3 + l)0'l + 32 + 3Z + 3)} 



-(Jl - 


+ -Ji 


- 1; 


) - 


(J2 - - 


"J2 


- 1) - (j 3 - -j3 - 1) 




+ -Ji 


- 1. 


J2 


-» -J2 


- r 


) + (j2 -> -J2 - 1,J3 -> -J3 


+0'i - 


+ -Ji 


- 1. 


J3 


-» -J3 


- r 


) 


"(ji - 


+ -Ji 


- 1. 


J2 


-> -J2 


- 1. 


J3 -> -J3 - !)• 



To obtain the sum of squares of eigenvalues with fixed j we subtract the 
result for f 2 (w, ji, J2, J3) with m = — j from the result for m = — j — 1 (it 
follows from the fact that the eigenvalues are the same for all m in (j, m) 
multiplets). We therefore expand (1 — w)f 2 (w,ji,j 2 ,j3) given in PU)) in w, 
and the coefficient in front of w - - 7 is the desired sum of squares of eigenvalues 
with fixed j. Hence for example: 

• the eigenvalue for j = (ji + j 2 + J3) is equal to 

• for j = (j! + j 2 + J 3 - 1) eigenvalues are equal to ±^Jjij 2 j 3 (ji + J2 + J3) 

• for j = (ji + j 2 + j 3 — 2) there is an eigenvalue and two eigenvalues 

equal to ±yj (4jij 2 j 3 - J1J2 - J2J3 - jij3)(ji + J2 + J3 - 1) + J1J2J3 (if 
any of the ji is equal to 1/2 then some of these eigenvalues may be 
missing). 

We may note that for the n-valent vertex [n > 4) the method works 
in exactly the same way: in the q 2 operator, beyond the already discussed 
terms, we encounter also terms of the form 

ji jj jk jl jm jn /91 \ 

^ijk°\° 2 "3 ^lmn \° 2 ^ ■ \' L ' L ) 

After sandwiching between |m 1; m 2 , m 3 , m 4 ) we have to assign k = 3, n = 3 
and therefore (i,j,k,l) = (+,—,—,+) or (i,j,k,l) = (—,+,+,—), so that 
the appropriate term in f 2 (w, ji, j 2 , J3, j'4) reads 

E ^ mi+m2+m3+m4 ^m 3 m 4 (ji(ji + 1) - ml) (j 2 (j2 + 1) - mf) 

mi,m2,m,3,m4 

(22) 

and the sum can be performed without difficulty. 

As can be easily seen the terms of the form eijkJ\J 2 J^e[ mn J\J™ J§ or 
tijkJ\J 2 Jz e imnJ\J™J§ give no contribution to f 2 . 
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We have also calculated along the same lines 

h(wJuj 2 Js) ■= Tr (q 4 w J * +J i +J l) (23) 
but the resulting formulae are too long to be presented here. 

5 The density distribution for eigenvalues 

It would be very desirable to have also higher-order moments, 

T r 9 \j,h,hj3J (24) 

since we could then reconstruct the full density by inverse Mellin transform 
of the relation 

'I'' <r" Jd\p 3 (X)X 2n . (25) 

Since q is an operator in a finite-dimensional space, the density distribution is 
just a sum of Dirac delta functions; we aim at finding a smooth approximation 
of this distribution for large values of the spins involved. 

In this section we will calculate higher moments (Tr g 4 , Tr q 6 and Tr q 8 ) 
as a function of j±, j'2, and j 3 , but with j not fixed. 

Introducing the notation 

(B) == If, (26) 



exp (-{A) + \{{A*) - {Af) 



--((A*)-3(A 2 )(A) + 2(A) 3 ) (27) 



we have the relation 

{e- A , - 

1 

3 

+^((A 4 ) - A(A 3 )(A) - 3(A 2 ) 2 + 12(A 2 )(A) 2 - 6(A) 4 ) + ... 

Neglecting terms of the third order and higher in the exponent (Gaussian 
approximation) we have 

(e-^> w e -<«9>+^ «? 2 )-«> 2 ) = | dAp G (A)e- aA (28) 
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and since (q) = we have 



where 

P = (?). (30) 

In the Gaussian approximation 

(q 2n ) = (2n-l)!!(? 2 ) n (31) 

As we already know, 



m = ^3, (32) 



where we introduced 



F = 1, 

Fl = Ji(ji + 1)+J 2 (j2 + 1)+J3(j3 + 1), 

F 2 = + l)j 2 (j2 + 1) + JlUl + l)j 3 (j3 + 1) + U]2 + l)j 3 (j3 + 1), 

F, = ./:(./. • M./21./2 • ■ I)- (33) 

Using the form (|TT|) we obtain (with the help of the MAPLE program) 

(<f) = J| (24F 3 - 8F 2 + llFx - 12F ) , (34) 

(<f > = tt^t (l440F 2 - 1440F 3 F 2 + 480F 2 2 + 2538F 3 F 1 - 1656^ 

-6909F 3 F + 1042F 2 + 5158F 2 F - 4815FiF + 3222F 2 ) (35) 



and 



(f) = f 4480F| - 8960F 2 F 2 + 8064F 3 F 2 2 + 19936^1^ 

w 7 99225 V 3 32 32 31 

-78032F 3 2 F - 2688F 2 3 - 38304F 3 F 2 Fi + 180204F 3 F 2 
+41853F 3 F 2 + 18480F 2 F! - 385701^^^0 - 91932F 2 F 
-36351F 2 F 2 + 691740F 3 F 2 + 329841F 2 FiF + 18081F 3 (36) 
-549288F 2 F 2 - 204552F 2 F + 521316F 1 F 2 - 285282F 3 ) . 



8 



The result for (q 2n ) has the form F 3 ■ P n -\ where P ra _i is a polynomial with 
alternating signs, homogeneous of degree (n — 1) in F 3 , F 2 , F 1 , F . 

We were unable to find a general expression for the coefficients of these 
polynomials, however, as will be shown in the next section, we were able to 
prove that for {q 2n ) the coefficient in front of F 3 , which is the most important 
term for large spins, is equal to ( 2 n+i)^2n+i)!! an< ^ f° r n > 1 it differs from the 
Gaussian approximation result (2n — l)!!(2/9) n . 



6 Classical distribution of the eigenvalues 

In this section we will find the classical distribution (both integrated p(A) and 
in the subspace of fixed total spin pj(A))of the eigenvalues of the operator q 
in the limit of large spins. In this limit we can disregard any commutators of 
J/ with itself (since they introduce lower power of spins) and the case boils 
down to the classical density. 

For the classical density we have, 

r"'"dA ft (A)A- = Wl5i» (JlXj2) . j3 f>. 
J-nmz J (47rJ° 

<5((J 1 + J 2 + J 3 ) 2 -j 2 ), (37) 

where we integrate over angular positions of Ji, J 2 and J 3 , and N is a 
normalization factor: 

N = j ^ ^ (( Ji + J 2 + J 3 ) 2 - J ) • (38) 

We introduce also the integrated density satisfying 

f jlhh a\ f\\ \2n f dQidQ 2 dn 3 2n 

/ dAp(A) A =/ — — — ((JixJ 2 )-J 3 ) . (39) 

J -313233 J (47TJ 

Since only relative angles are important, we can fix some positions: 

fl 3 = 0, 2 = (40) 

and then for example 

( Ji x J 2 ) • J 3 = j 1 j 2 j 3 sin 6>i sin 9 2 sin <p. (41) 
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We introduce coordinates on S 3 : 



x = sin 9 1 sin 9 2 sin 0, 

y = sin 6\ sin 9 2 cos 0, 

z = sin #1 cos #2, 

t = cos^i. (42) 



Then 

(Jij 2 j3x) 2n 5 (jI + j 2 2 + j 3 2 + 2 3l32 t + 2j 1 j, 



+2J2J3 —2= + 2j 2 j 3 y - J 2 J • (43) 

Comparing the sides we obtain the final result, i.e. the classical distribu- 
tion function of eigenvalues in terms of j±, j 2 , j 3 , j: 

p,(A) = —i /^^25(A 2 + , 2 + , 2 + t 2 -l) 

5 (]l + j! + j! + 2 Jlj2 t + 2^ -^== 

zt \ 
+2j 2 j3^== + 2j 2 j 3 y - fj , (44) 

where 

A A (./I./2./V,!. (45) 

We can integrate over t: 

1 /-v 7 ^ /-2- r drd0 

Pi (A) : 



87rN 3l32 j 3 J Jo v / r2 + T 2v / 1 _ r2 _p 



E 5 I ii + jl + jl + Ijihe VI - r 2 - A 2 + 2j lj3 

c=±l 



r sin 



. r sin <W 1 — r 2 — A 2 

+2j 2 j3e , + 2j 2 j3^ COS0 — j' 

Vr 2 + A 2 



y/r 2 + \ 2 
\ 

(46) 
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It is difficult to perform the integral in this form. However, performing the 
same steps for the integrated distribution function we obtain a surprisingly 
simple result: 



1 rV i-A 2 rc l r 1 
p(X) = / - = arccos (A) (47) 

2jlJ 2 j3 -A) J r 2 + _ r 2 _ ^2 2j!j 2 j3 



for |A| < 1 and for |A| > I. 

Therefore we get the coefficient in front of Fg for the integrated trace 

(q 2n ) -> / dx- arccos(x)x 2n F" = y — p -K 1 . (48) 

w 1 J-i 2 v ' 3 (2n + l)(2n + l)!! 3 

Hence for example we recover the leading coefficients of (|32|). (|34|) . and 

tf.)^,, <«^F|, <4^||* (49) 

Concluding, by using the methods described in this paper we can learn 
something about the distribution of the eigenvalues of the volume operator in 
the quantum case, and there exists an explicit formula for this distribution in 
the limiting case of large spins, i.e. the classical case. It would be interesting 
to compare these results with the numerical analysis for low spins (up to 50) 
given in [5]. 

Besides being mathematically interesting, these results can be important 
for example in the statistical approach to loop quantum gravity. 

Several conceptual problems arise in this context by the results of the 
present paper. The first one is what observables (besides the volume) should 
be specified to actually describe a macroscopic state in a given region. Since 
the contribution to the volume from a given 4-valent vertex grows like j\ -ji-jz 
it is important to know whether these observables provide suppression of 
the contribution from very large spins rendering final expectation value of 
the volume finite. The second problem is connected with the first one and 
concerns the bulk entropy: there are many sets of graphs with many different 
assignments of spins on legs that give (macroscopically) the same volume. If 
we specify both the volume and the other observables we are left with a 
certain number of graphs. The question is, whether the logarithm of this 
number is connected with the bulk entropy much as sequences of spins are 
connected with the black hole surface entropy jU]. 
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